The R 2 phase-diagram of QEG and its spectral dimension 
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Within the gravitational asymptotic safety program, the RG flow of the R 2 truncation in three 
and four spacetime dimensions is analyzed in detail. In particular, we construct RG trajectories 
which emanate from the non- Gaussian UV fixed point and possess long classical regimes where 
the effective average action is well approximated by the classical Einstein-Hilbert action. As an 
application we study the spectral dimension of the effective QEG spacetimes resulting from these 
trajectories, establishing that the picture of a multi-fractal spacetime is robust under the extension 
of the truncated theory space. We demonstrate that regimes of constant spectral dimensions can 
either be attributed to universal features of RG fixed points or singular loci in the j3 functions. 
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1. Introduction 

Weinbergs asymptotic safety scenario [1, 2] strives 
toward constructing a quantum theory of gravity 
based on a non-Gaussian fixed point (NGFP) of the 
gravitational renormalization group (RG) flow (see 
[3, 4, 5] for reviews). This NGFP is supposed to 
control the behavior of the theory at high energies 
and renders it safe from unphysical divergences. In 
order to ensure the predictivity of the construction, 
the NGFP has to come with a finite number of UV- 
relevant directions, i.e., the space of RG trajectories 
attracted to the fixed point at high energies is finite- 
dimensional. Moreover one has the phenomenologi- 
cal requirement that at least one of the RG trajecto- 
ries emanating from the NGFP possesses a "classical 
limit" where the Einstein-Hilbert action constitutes 
a good approximation of the gravitational interac- 
tion. The resulting quantum theory for gravity is 
called Quantum Einstein Gravity (QEG). 

A central element for investigating Asymptotic 
Safety is the functional renormalization group equa- 



tion (FRGE) for the gravitational average action 
[6] , which allows the systematic construction of non- 
perturbative approximations of the theories ft func- 
tions. Projecting the RG flow to subspaces of suc- 
cessively increasing complexity, the existence of a 
NGFP has been established in the Einstein-Hilbert 
truncation [7, 8, 9], the R 2 truncation [10, 11], f(R) 
truncations [12], and truncations including a Weyl- 
squared term [13]. Moreover, refs. [9, 14] studied 
the properties of the NGFP in spacetime dimension 
more than four, while the quantum effects in the 
ghost sector have been investigated in [15]. Further- 
more boundary terms are considered in [16]. QEG 
with Lorentzian signature was investigated in [17] 
and an computer based algorithm for evaluating the 
flow equations was proposed in [18]. Besides provid- 
ing striking evidence for the existence of the NGFP 
underlying Asymptotic Safety, the works including 
higher derivative interactions also provided strong 
hints that the fixed point possesses a finite number 
of relevant directions. 

While there is solid evidence in favor of a NGFP 
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Figure 1: Classification of the RG flow in the 
Einstein-Hilbcrt truncation [20]. 



in quantum gravity, a lot less is known about the 
IR physics of the RG trajectories emanating from 
it. In case of the Einstein-Hilbcrt truncation the 
RG trajectories have been classified in [9] (also see 
[19] for a recent account) and are shown in Fig. 1. 
All trajectories with positive Newtons constant are 
attracted to the NGFP in the UV and develop a 
classical regime close to the A axis. Based on this 
classification, the RG trajectory realized by Nature 
has been identified in [21]. 

Since the [3 functions of more sophisticated trun- 
cations quickly become very involved, there is only 
very little known about its dynamics of the RG flow 
beyond the Einstein-Hilbert case, however. One aim 
of the present work is to fill this gap by constructing 
families of trajectories originating from the full flow 
equations of the R 2 truncation in d = 3 and d = 4 
spacetime dimensions. 

Our motivation for this study is two-fold. With re- 
spect to the fundamental aspects underlying Asymp- 
totic Safety, it is important to understand to which 
extend the features of the RG flow shown in Fig. 1 
remain robust when the effect of additional running 
couplings are included. In particular the existence of 
a classical regime can only be established by study- 
ing the dynamics of the RG flow and the interplay of 
different fixed points. Surprisingly, our analysis of 
the R 2 truncation also provides further insights on 
the IR fixed point recently discussed in [22, 23, 19]. 

A more phenomenological motivation arises from 
studying quantum effects in gravitational systems 
via an "RG improvement", taking into account the 



scale-dependence of the gravitational coupling con- 
stants when studying, e.g., black holes [24, 25], cos- 
mology [26, 27, 28, 29, 30, 31, 32, 33, 34], or galaxy 
rotation curves [35]. The RG improvement requires 
rather explicit knowledge of the RG trajectories un- 
derlying the physical system. Thus our results may 
be used to investigate the robustness of previous 
studies. 

As a first application in this direction, we ex- 
tend the formalism [36, 37] to study the multi- 
fractal properties of the effective QEG spacetimes 
captured by the spectral dimension D s . The in- 
teresting question from the functional renormaliza- 
tion group point of view is whether the fractal prop- 
erties found within the Einstein-Hilbcrt truncation 
[37] stay valid in higher truncations and if there are 
new features special to the higher derivative ansatz. 
In this sense, the R 2 truncation constitutes a natural 
choice for computing on shell corrections to the frac- 
tal properties of spacetime observed in the Einstein- 
Hilbert case. 

Notably, the interest in spectral dimension of 
spacetime is not limited to the functional renormal- 
ization group, since it is also accessible within other 
approaches to quantum gravity, as e.g., causal dy- 
namical triangulation (CDT) [38], euclidean dynam- 
ical triangulation (EDT) [39] , loop quantum gravity 
and spin foam models [40] or in the presence of a 
minimal length scale [41]. Studying the spectral di- 
mension from various different angles may unravel 
universal features in the quantum structure of space- 
time. On more phenomenological grounds one can 
build toy models of spacetime which encompass spe- 
cial features of the spectral dimension [42] . Among 
further developments a fractional differential calcu- 
lus [43] was used in [44, 45] to assemble spacetimes 
with the fractal features observed in quantum grav- 
ity. Thus D s constitutes a quite useful probe for 
quantum gravity effects. 

The remaining parts of the paper are organized 
as follows. The phase diagrams of the R 2 truncation 
are constructed in section 2. Section 3 introduces the 
formalism for investigating the spectral dimension 
D s and we give explicit examples in section 4. We 
close with a discussion of our results in section 5. 
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2. Fixed points, phase diagrams and Asymptotic 
Safety 



When studying the RG flow arising from (2.3), we 
first note 



The gravitational part of the R 2 truncation takes 
the form 



pgrav 
k ' 



d d Xy /g 



(-R + 2\ k )+±R 2 



(2.1) 



The ft functions governing the dependence of the di- 
mensionful Newtons constant G k , cosmological con- 
stant Afe and R 2 coupling b k on the RG scale k have 
been derived in [11] for general dimension d and 
regulator lZ k . In order to facilitate our numerical 
analysis, we will work with the optimized cutoff [46] 
throughout. 

2.1. /? functions of the R 2 truncation 

In order to study the RG flow resulting from the 
ansatz (2.1), it is convenient to work with the di- 
mensionless couplings 



9k 



G k 



A,, 



bk 



(2.2) 



Their scale dependence is governed by the /3 func- 
tions 



d t gk = P g (g,\,b) , 
d t Xk =0\{g,\,b), 
d t b k = Pb(g,\b), 



(2.3) 



with t = In k being the RG time. Their explicit form 
can be found in [11] to which we refer for further 
details. 

There the (3 functions (2.3) have been derived with 
the help of the transverse-traceless decomposition of 
the fluctuation fields. This decomposition give rise 
to "zero-mode" contributions proportional to <5<2.4- 
A detailed investigation shows that these terms give 
rise to additional singularities of the /3 functions, 
which also have a significant influence on the prop- 
erties of the RG flow in the UV (see below for a 
more detailed discussion). Since there is no reason 
to expect that similar singularities arise when the 
computation is carried out without the use of the 
transverse traceless decomposition, we remove these 
contributions in four dimensional calculations. 1 



1 Alternatively, this can be seen as analyzing the "regular- 
ized" system of ft functions (2.3) in d = 4 — e dimensions. 



/3 ff (A,<7,&)| s =o=0, &(A,<7,&)| 6 =o=0, (2.4) 

which implies that RG trajectories cannot cross the 
,9 = plane and the b = plane. As a consequence, a 
RG trajectory starting with a positive Newtons con- 
stant cannot pass to a regime with negative Newtons 
constant dynamically. Eq. (2.4) does not ensure the 
positivity of bk, however, since the coupling b k can 
change sign by passing 1/bk — 0, which is a regular 
point of the (3 functions. 

In addition, the j3 functions (2.3) contain a singu- 
lar locus where the denominators of the anomalous 
dimensions tjn = dt lng/c and r\h = dt ln&fc vanishes. 
For d = 4, the resulting singular locus is depicted in 
Fig. 2. 

The final ingredient in understanding the RG flow 
of the theory are the fixed points where (3 a (u*) = 
for all a = 1, 2, • • • . In the vicinity of such a fixed 
point of the RG equations (2.3), the linearized flow 
is governed by the Jacobi matrix B = (B a7 ), B ai = 




Figure 2: The singular locus in theory space where 
the anomalous dimensions t/n and 775 
diverge. 
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d 7 (3 a (u*): 

kd k u a (k) = y~] B a7 (uj(k) - u*) . (2.5) 

7 

The general solution to this equation reads 

«„(*)=<+ 2 i^y , (2.6) 

where the V 1 's are the right-eigenvectors of B with 
eigenvalues —0/, i.e., ^ 7 -B a7 V^ = — 0/V^. Since 
B is not symmetric in general the 0/'s are not guar- 
anteed to be real. 

2.2. The phase diagram in d = 4 

In d = 4, and for an exponential cutoff the fixed 
point structure of the system (2.3) has already been 
analyzed in [11]. Besides the Gaussian fixed point 
(GFP) one encounters a NGFP at {X*,g*,b*} = 
{0.330, 0.292, 183.5} with three UV-relcvant eigendi- 
rections 

0i,2 = 2.15 ±3.79i, 03 = 28.8. (2.7) 

In order to analyze the corresponding RG flow on 
phase space, we first repeat this analysis using the 
optimized cutoff. This also gives rise to a GFP and 
a NGFP. 

Gaussian Fixed Point 

Firstly, the system (2.3) exhibits a GFP, 

GFP: {X*,g*,b*} = {0, 0, 0}. (2.8) 

Its stability coefficients are 

{01,02,03} = {2,0,-2}, (2.9) 

and thus correspond to the mass dimensions of the 
couplings. The linearized stability analysis shows 
that we have one IR-attractive, one IR-repulsive and 
one marginal eigendirection, in agreement with stan- 
dard power counting. Since the GFP lies on the sin- 
gular locus shown in Fig. 2, the corresponding stabil- 
ity matrix and its eigendirections depend on the pre- 
cise order in which the limit {X*,g*,b*} -> {0,0,0} 
is taken. To get a rough idea of the behavior close 
to the GFP one can build the following limits inde- 
pendent of the approaching direction 

<9tAfc| s= o,&=o = -2\ k + (D(\l), 



d t g k \\=o,b=o = 2g k + 0(g 2 ), (2.10) 

Since d t bu is quadratic in b it is IR-repulsive for pos- 
itive b and IR-attractive for negative b. 
Non-Gaussian Fixed Point 

In addition the four-dimensional RG flow gives 
rise to a unique NGFP, which is located at 
{A* ,#*,&*} = {0.170,0.754,336.2}. Its stability 
coefficients are 

0i= 3.95 , 02,3 = 1-56 ±3.31i. (2.11) 

The positive real part of the critical exponents indi- 
cates that any RG trajectory in the vicinity of the 
fixed point is attracted to the NGFP in the UV while 
the non- vanishing imaginary part of 02,3 shows that 
the trajectories will spiral into the NGFP. 

Notably, the critical exponent 9\ is significantly 
smaller than the one found in eq. (2.7). This dif- 
ference can be traced back to the zero-mode con- 
tribution in the (3 functions. Including these the 
optimized cutoff actually gives rise to two NGFPs 

NGFPi : {A* , g* , b* } = {0.163, 0.744, 345} 

NGFP 2 : {A* , g* , b* } = {0.173, 0.696, 455} 

(2.12) 

with 

NGFPi : 0i = 3.49 , 2 3 = 1-65 ± 3.RH , 

(2.13) 

NGFP 2 : 0i = 25.6, 2 , 3 = 1-40 ± 2.78 . 

The critical exponents of the second one are very 
close to the one found in the analysis [11]. When re- 
moving the zero-mode contribution NGFP 2 vanishes 
and the NGFPi falls on top of the NGFP (2.11). Re- 
markably it was found in [12] that the largest criti- 
cal exponent at the NGFP is of order 25 for the R 2 
truncation and approximately 2 when the R 3 term 
is taken into account as well. Thus our FP with the 
largest critical exponent being 0i = 3.95 fits well to 
the higher order results. 
Trajectories 

Starting at the NGFP there are several types of tra- 
jectories flowing toward the IR. In Fig. 3 we show a 
particular sample of them which exhibit a crossover 
to the GFP and are thus likely to give rise to a 
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Figure 3: A sample of possible trajectories including 
the examples A and B. The fixed points 
are marked with dots. 



classical regime. In analogy to the phase diagram 
of the Einstein-Hilbert truncation depicted in Fig. 
1, we distinguish RG trajectories of Type la and 
Ilia, according to their IR behavior as follows: Upon 
their crossover, trajectories of Type la run toward 
negative values A. For k — > 0, they approach the 
point {A, g,b} — {— oo,0,0}. They exhibit a classi- 
cal regime with a negative value of the dimensionful 
cosmological constant X k . In contrast, the trajecto- 
ries of Type Ilia leave the GFP regime in direction of 
positive cosmological constant. They may first flow 
parallel to the A axis toward increasingly positive 
values Afc. Here the trajectories develop a classi- 
cal regime with a positive cosmological constant X k . 
Before hitting the singular locus depicted in Fig. 2, 
however, the trajectories essentially turn perpendic- 
ular to the A-g-plane and flow toward increasing val- 
ues bk- This turn can already occur at rather small 
values Afe, such that the flow is along the b axis. 
Trading the coupling bk for 1/bk, shows that the R 2 
coupling changes sign. At the end all these trajec- 
tories tend toward the point {A, g, b} = {1/2, 0, 0~}. 
Note that these observations are in perfect agree- 
ment with the recent suggestions [22, 23, 19] that 
the corresponding point {A, g} — {1/2,0} in the 
Einstein-Hilbert truncation actually constitutes an 



IR fixed point. 

For later purposes, we marked two explicit Type 
Ilia sample trajectories A and B. These arise from 
the initial conditions 

A : {A ini , g ini , & ini } = {0.4999 , 10" 10 , 10 9 } , 

B : {A ini , flw , &ini } = {1(T 5 , 10- 8 , 500} . 

(2.14) 

They will form the basis for our discussion of the 
spectral dimension in section 4. 

Finally, we remark that we have not been able to 
construct a RG trajectory that connects the NGFP 
with the GFP and which would constitute the ana- 
logue of the separatrix (Type Ha trajectory) shown 
in Fig. 1. The non-existence of a separatrix in the R 2 
truncation is caused by the singularity structure of 
the (3 functions: The RG flow away from the NGFP 
cannot be aligned to the IR-attractive eigendirection 
of the GFP. Thus we can finetune the flow such that 
the trajectory passes arbitrarily close to the GFP, 
but there is no solution connecting the fixed points. 

2.3. The phase diagram in d = 3 

In this subsection we will repeat the discussion of the 
previous one for the case of three spacetime dimen- 
sions. The three dimensional case is a model very 
well suited for comparison to Monte Carlo data as 
e.g. [47] since the computations in three dimensions 
are less expensive within the CDT framework. 
Gaussian Fixed Point 

Again, the GFP is situated at {A*, #*,&*} = 
{0,0,0}. Its associated stability coefficients are 
given by 

{01,02,0a} = {2,1,-1}. (2.15) 

These correspond to the mass-dimension of the (di- 
mensionful) coupling constants. Eq. (2.6) then in- 
dicates that there are two IR-repulsive and one IR- 
attractive eigendirection. These depend on the di- 
rection of the approach to the GFP which is again 
situated on the singular locus. The limits 

<9tAfc| g= o,b=o = — 2Afc + 0(X k ) , 
d t g k \ x = ,b=a = 9k + 0(g 2 k ), (2.16) 
d t b k \ g=0 ^x=o = — bk — ■p^pgffrfc ' 
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again are independent of the direction of the ap- 
proach and give a rough idea of the behavior close 
to the GFP. This time the 6 direction is IR-repulsive 
for positive and negative 6 values. 
Non- Gaussian Fixed Points 
Investigating the limit g = 0, one first encounters 
a non-trivial FP at {A*, g*, b*} = {0, 0, -f|(47r) 2 }. 
This FP has two eigendirections which both have 
no g component. Thus all trajectories connected to 
this FP can not leave the g — plane. We do not 
consider this FP any further. 

In addition the system (2.3) gives rise to a 
physical NGFP which is located at {X*,g*,b*} = 
{0.019, 0.188, 126.4}. Its critical exponents read 

{0%, O2, 3 } = {8.39, 1.86, 1.35} . (2.17) 

Thus this NGFP comes with three UV attractive 
directions and all trajectories in its vicinity are 
dragged into this FP in the UV. The reality of the 
critical exponents indicates that the curling into the 
NGFP seen in four dimensions is absent. By vary- 
ing the dimension in our theory one can see that the 
NGFPs (2.17) and (2.11) are continuously connected 
in d. 

Besides this NGFP one finds two additional un- 
physical fixed points (UNGFPs). For completeness, 
their position and stability properties are collected 
in Table 1. By tracing the flow of the fixed points 




Figure 4: Position of the UNGFPi (blue) and 
UNGFP 2 (red) depending on d. The two 
fixed points annihilate at d = 3.84. 



UNGFPi and UNGFP 2 under an increase of the di- 
mension d, one finds that they annihilate each other 
for d = 3.84 at {X*,g*,b*} = {0.27,0.54,76.9} as 
shown in Fig. 4. A similar continuous treatment of 
the dimension for the Z2-effective potential in [48] 
showed the appearance of new fixed points when de- 
creasing the dimension d as well. 

Compared to the NGFP the UNGFPs are situ- 
ated on the other side of the singularity locus. In 
principle, both UNGFPs do act as UV attractors 
for the RG trajectories spanning their UV-critical 
hypersurface. Tracing the RG flow toward the IR 
reveals, however, that the corresponding RG trajec- 
tories do not flow to a region in theory space, where 
General Relativity provides a good approximation. 
Thus, while still giving rise to RG trajectories which 
are asymptotically safe in the UV, they do not give 
rise to a classical regime. 
Trajectories 

In three dimensions the phase diagram is similar to 
the one in four dimensions. Again we focus on RG 
trajectories that exhibit a crossover from the NGFP 
to the GFP. Some typical examples are depicted in 
Fig. 5 and for later purposes we marked the trajec- 
tories A and B. These are obtained with the initial 
conditions 

A : {A init , ffinit , &init } = {0.49 , 10~ 8 , 5 x 10 4 } , 
B : {A init , g init , 6 init } = {10- 5 , 10- 8 , 500} . 

(2.18) 

Besides the GFP and the NGFP, Fig. 5 also de- 
picts the two UNGFPs. These are separated from 
the NGFP by the singular locus shown in Fig. 2. 
Moreover, we highlighted the prospective IR fixed 
point {X,g,b} = {0.5, 0,0~}. 

The first obvious difference to the four dimen- 
sional plot in Fig. 3 is that in three dimensions the 
approach toward the NGFP is straight since the crit- 
ical exponents (2.17) are real, compared to the com- 
plex ones (2.11) in four dimensions. 

As we did in four dimensions we can divide the 
trajectories into Type la and Type Ilia trajectories. 
The former are going toward {A, g, b} = {—00, 0, 0} 
in the IR. The latter pass the GFP and head toward 
infinite b values at positive Xk afterward. Again the 
use of 1/6 instead of b shows that the sign changes 



(i 





A* g* b* 


61 6 2 #3 


UNGFPi 
UNGFP2 


0.364 0.147 18.167 
0.099 0.216 18.835 


1.56±4.84i -9.85 
0.19±0.97i 1.68 



Table 1: Additional fixed points of the beta functions (2.3) in d = 3. These fixed points are labeled "unphysical" 
since the RG flow emanating from them does not give rise to a classical regime where General Relativity 
constitutes a good approximation. 



and the Type Ilia trajectories go toward the point 
{A, g, b} = {1/2, 0,0~} with negative b values. 

In three (as well as four) dimensions we get large 
contributions to the j3 functions between the GFP 
and the NGFP as well as close to A = 1/2. These are 
caused by the singularity structure of the truncated 
theory space and lead to a strong running in the b 
direction. 



0.25 




Figure 5: Phase diagram in three dimensions with 
the GFP, the NGFP and the two 
UNGFPs. We depict a set of trajecto- 
ries which develops a long classical regime 
and also includes the examples (2.18) high- 
lighted with bold black lines. 



3. RG-improved diffusion processes 

As a first application, we use the phase-space study 
of the last section to investigate the influence of 
higher-derivative terms on the fractal-like properties 
of the effective QEG spacetimes. In this course, we 
start by adapting the "RG improvement" scheme for 
the diffusion of a test particle on the QEG spacetime 
[36, 37] to the higher- derivative action (2.1). 

The key idea of the RG improvement is to study 
diffusion processes with a spacetime metric, which 
depends on the RG scale k 



d T K(x, x'\ T) = -A(k)K(x, x'; T) . 



(3-1) 



Resorting to the "flat-space" approximation, the 
scale k is identified with the momentum p of the 
plain waves used to probe the structures of space- 
time and thus provides a (inverse) length scale at 
which the system is probed. Taking the trace of 
K(x,x';T), gives the average return probability 



P(T) 



d d k 
(2^p 



-k 2 F(k 2 )T 



(3.2) 



where the function F(k 2 ) relates the Laplacian at 
a fixed reference scale ko (taken in the IR) and the 
scale k: A(k) = F(k 2 )A(ko). The scale-dependent 
spectral dimension is defined as the logarithmic T 
derivative 



D s (T),-2« 
V ' d\nT 



(3.3) 



For F(k 2 ) = 1 this relation reproduces the classical 
case D s (T) = d; thus all quantum corrections to the 
spectral dimension are encoded in F(k 2 ). 

Assuming that F(k 2 ) oc k s undergoes power law 
scaling, the integral (3.2) can be evaluated analyti- 
cally, yielding P(T) oc T'W+V. Substituting this 
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relation into (3.3) yields the spectral dimension 

2d 



This equation allows to read off the function F{k 2 ) 
relating the metrics at the scale k and ko 



V S {T) 



(3.4) 



2 + 6 ' 

For the Einstein-Hilbert case [37] 

<5 EH (. 9 ,A) = 2 + A- 1 /3 A (. 9 ,A) (3.5) 

can be expressed as a function on theory space. The 
main task of this section is to derive the analogue of 
this expression for the ansatz (2.1), taking the R 2 
contribution into account. 

Following [37], we first write down the equations 
of motion arising from (2.1). Upon integrating by 
parts and dropping the surface terms these are 

+ 2X k + 16 ^ k kR2 ) 9^ + ( 2 - 64w £ kR ) Rfl " 
+ ~^ (D^D U R - (D 2 R) g^ v ) = . (3.6) 

In contrast to the classical equations of motion, the k 
dependence of the coupling constants promotes (3.6) 
to a one-parameter family of equations of motion, 
each yielding an effective description of the physics 
at the fixed scale k. In order to extract the k depen- 
dence of the metric, we first keep the scale k fixed 
and substitute the ansatz 



<r\ k = F(k)<r\ k0 



(3.11) 



with F(k) = X k /X ko . In a regime where F(k 2 under- 
goes power-law scaling for an extended fc-interval, 
the function 6(k) can then be obtained as 5(k) = 
kdk In F{k) and reads 



%, A, 6)| d = 



2 + A" 1 / 



(3.12) 



Formally, this result has the same form as in the 
Einstein-Hilbert case (3.5). For the R 2 truncation 
the /? function (i\ depends on the three couplings 
A, g, b, however, so that there are non-trivial correc- 
tions. 

For d ^ 4 the situation is slightly more compli- 
cated. In this case the solutions of the quadratic 
equation are 



± (d-2)b k 
l 'k 327r(d-4)G fc 



where 



h d 



128(i(d-4)7r 



k Afc 



(3.13) 



(3.14) 



R^v {g\k) — -j9tj.u\k ■ 



(3.7) 



which implies constant curvature R(g\k) = c k . For 
this ansatz, the second line in (3.6) is zero while the 
first line yields 

2X k -\{d-2)c k + \ (d-4) i^ c 2 fe = 0. (3.8) 

At this stage, it is convenient to distinguish the two 
cases d = 4 and d ^ 4. 

For d = 4 the term quadratic in c k vanishes and 
cq. (3.8) is easily solved 



(d-2) 2 ■ 

Investigating the conditions implied by eq. (3.10), 
one then observes that may not be positive def- 
inite. In particular for d = 3 and positive G and b, 
c k is negative. Therefore we discard this solution 
and focus on . Notably, also has a well-defined 
limit <i — ^4 where it reduces to (3.9). Thus select- 
ing the CjT branch will lead to a spectral dimension 
which is continuous in d. 

The function F(k 2 ) originating from cjT can then 
be obtained completely analogous to the case d = 4 
and reads 



F(k) = 



Gob k 1 - yjl - h d G k \ k /b k 



(3.15) 



c k \d=4 = 4Afc . 



(3.9) 



Comparing solutions (3.7) at two different scales k 
and ko, where ko is a fixed reference scale taken in 
the classical regime, and using the identity R£(cg) = 
c~~ 1 R£(g), which holds for constant c > 0, then 
yields 

R;(g k ) = R;(^g ko ) . (3.10) 



Gkb Q l - y/i - h d GoXo/b ' 

Again, this leads to a S(k) which is given by a func- 
tion on theory space 



9 b (1 - VI - h d g\/b) y/1 - h d g\/b 

(3.16) 

Substituted into the spectral dimension (3.4), this 
result captured the R 2 corrections to D s . Its prop- 
erties will be investigated in the next section. Keep 
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in mind, however, that D s is ill-defined if 5 changes 
rapidly [37]. 

4. Fractal properties of QEG spacetimes 

In this section, we analyze the spectral dimension 
obtained along the sample trajectories (2.14) and 
(2.18) in d = 4 and d = 3, respectively. For this 
purpose we evaluate (3.16) along the trajectory and 
substitute the result into (3.4) to obtain the spectral 
dimension D s (t) depending on the logarithmic RG 
time t — In k. 

Spectral Dimension in d = 4 

For d = 4 the two exemplary trajectories A and B 
are depicted in Fig. 3 . Both trajectories emanate 
from the NGFP in the UV and flow toward the GFP. 
After passing the GFP trajectory A flows toward 
A = 1/2 along the A axis before running toward large 
values bk- Trajectory B leaves the GFP regime close 
to the b axis instead. 

The spectral dimension along trajectory A is 
shown in Fig. 6. The diagram displays three dis- 
tinguished plateaus marked by the solid blue lines. 
From the UV (large t) to the IR (small t) these are 
the NGFP plateau with D s = 2, the semi-classical 
plateau with D s ~ 1.5 and the classical plateau 
where D s — 4, respectively. They are connected by 
interpolating regimes (dashed), where 5(t) changes 
rapidly and thus D s (t) is not well defined. The poles 
in these regions are caused by S(i) — > —2, which in 
turn is linked to A& changing sign when passing from 



the NGFP to the GFP. 

Surprisingly, the spectral dimension of the semi- 
classical plateau, D s — 1.5 does not agree with 

the one observed in the Einstein-Hilbert truncation 
£>eh = 4/3 j 37 ] Thig difference 

can be traced back 
to the fact, that in the R 2 truncation the flow for this 
intermediate part of the trajectory is dominated by 
the singular locus Fig. 2. This locus induces a scaling 
behavior, k s , S ~ 3.3, different from the fc 4 -scaling 
associated to the GFP. 

The spectral dimension of trajectory B, shown 
in Fig. 7 also exhibits three plateaus which are 
connected by short interpolating pieces where 5(t) 
changes rapidly. The NGFP plateau (D s = 2) and 
the classical plateau (D s = 4) have the same spec- 
tral dimension as for the trajectory A. In this case, 
the semi-classical plateau forms at D s = 4/3, how- 
ever, and thus has the same value as in the Einstein- 
Hilbert truncation. In terms of the underlying RG 
flow it is created by the trajectory flowing away from 
the GFP along the b axis (see the Appendix for an 
analytic investigation). Notably, there is no plateau 
corresponding to the scaling caused by the singular 
locus. This feature is shadowed by the poles where 
5 ->■ -2. 

Spectral Dimension in d = 3 

When investigating the spectral dimension in d = 3, 
we first notice that the functional form of (3.16) 
is quite different from the four-dimensional case. 
In particular S(X, g, b) now depends on the /3 func- 
tions of all three couplings. In order to see if the 
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Figure 6: Spectral dimension D s depending on RG Figure 7: Spectral dimension D s depending on RG 
time t for the trajectory A in d — 4. time t for the trajectory B in d = 4. 
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new structures still support a multi-fractal picture 
of spacetime, we again analyze D s (t) along the two 
exemplary trajectories A and B given in Fig. 5. 

The spectral dimension resulting from trajectory 
A is depicted in Fig. 8. Also here we observe three 
plateaus (solid line) which are connected by short 
interpolating pieces (dashed line) where the spec- 
tral dimension is not well-defined. In the IR and in 
the UV D s (t) again has the expected behavior: at 
small RG times t it equals the spacetime dimension 
three while for large RG times the NGFP induces 
a spectral dimension which is half of the spacetime 
dimension, D s = 3/2. 

The semi-classical plateau has D s ~ 1.3 and again 
develops during the crossover between the NGFP 
and the GFP. The poles framing the plateau are 
caused by A& changing sign which again leads to 6 = 
—2. Similarly to the four-dimensional case, the value 
of the semi-classical plateau does not correspond to 
the one observed in the Einstein-Hilbert case Df" = 
1. Instead of a property of the GFP this plateau is 
created by scaling induced by the singular locus and 
thus has a quite different origin. 

A special feature occurs when studying the spec- 
tral dimension along trajectory B, which is shown 
in Fig. 9. Besides the classical and NGFP plateaus, 
this trajectory develops two semi-classical regimes: 
The plateau with D s ~ 1.3 is caused by the trajec- 
tory passing close to the singular locus when flow- 
ing from the NGFP to the GFP. The second plateau 
with D s ~ 0.85 appears during the flow away from 
the GFP close to the b axis toward the classical 
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Figure 8: Spectral dimension D s depending on RG 
time t for the trajectory A in d = 3. 



regime. The latter plateau is the one that, in d = 4, 
corresponds to the semi-classical plateau found in 
the Einstein-Hilbert truncation. Thus it is natural 
to identify the shift Df H =14 Df 2 ~ 0.85 as the 
correction induced by the inclusion of the i? 2 -terms 
in the computation. 

The four different scaling regimes underlying Fig. 

9 are also visible in the running of the absolute value 
of the dimensionful cosmological constant A in Fig. 

10 whereas the running of the dimensionful New- 
ton constant shows only two distinct scaling regimes. 
Moreover the running of A shows two peaks indicat- 
ing the change of the sign. 

5. Conclusion 

In this work, we used the j3 functions [11] to con- 
struct the phase diagram of Quantum Einstein Grav- 
ity (QEG) in the R 2 truncation. For d — 4 and d = 3 
the resulting RG trajectories are shown in Figs. 3 
and 5, respectively. Notably, the flow is governed by 
the interplay of its fixed points. In this sense the 
phase diagram of the R 2 truncation is quite sim- 
ilar to the one of the Einstein-Hilbert truncation: 
the trajectories with a long classical regime result 
from a crossover from the non-Gaussian fixed point 
(NGFP) controlling the UV behavior of the theory 
to the Gaussian fixed point (GFP) responsible for 
its classical properties. Remarkably, the value of the 
R 2 coupling is easily compatible with experimental 
bounds [49, 50]. 

As an important new feature, our phase diagram 
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Figure 9: Spectral dimension D s depending on RG 
time t for the trajectory B in d — 3. 
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Figure 10: Absolute value of dimensionful cosmolog- 
ical constant A and dimensionful Newton 
constant G depending on logarithmic RG 
scale t along the exemplary trajectory B. 



is molded not only by its fixed points, but also by 
singular loci where the anomalous dimensions di- 
verge. The two regimes where this new mechanism 
is at work are the crossover between the NGFP and 
GFP and the flow close to A = 1/2. In both cases 
the singularities induce a strong running of the R 2 - 
coupling b. Especially close to the A = 1/2 locus the 
RG flow is essential parallel to the A = 1/2 boundary 
and the strong b running allows to avoid this singu- 
larity. The resulting flow pattern closely resembles 
the one expected from the IR fixed point recently 
discussed in [22, 23, 19]. 

As a first application of our classification, we stud- 
ied the spectral dimension D s of the effective QEG 
spacetimes. This is the first instance where higher- 
derivative corrections are taken into account in an 
RG- improvement scheme. As a central result, we 
establish that the multi-fractal structure of the ef- 
fective spacetimes [37], reflected by several plateau 
values for D s , is robust upon including R 2 effects. 
All our examples give D s = d in the classical regime 



and D s = d/2 in the UV. The latter constitutes the 
universal signature of the NGFP in the spectral di- 
mension [36]. 

Between the classical and UV regime we obtain 
the "semi-classical" plateau which, in the Einstcin- 
Hilbert case, is closely linked to universal properties 
of the RG flow close to the GFP and situated at 
D s = 4/3 (d = 4) and D s = 1 (d = 3). The R 2 
construction reveals a second mechanism that can 
lead to a similar "semi-classical" plateau: scaling 
induced by singularities. For our examples this in- 
duces semi-classical plateaus with D s ps 3/2 (d = 4) 
and D s ps 1.3 (d = 3). Thus the semi-classical 
regime is sensitive to the details of the underlying 
RG trajectory. In this light, it is not surprising 
that different Monte-Carlo studies of the gravita- 
tional path integral report different plateau values 
for D s at intermediate scales [47, 39]. 
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A. The GFP in d = 4 

In this section we investigate the GFP in four space- 
time dimensions a little bit further. Since we are 
interested in the behavior of the trajectory B in Fig. 
2 we evaluate the stability matrix by first taking the 
limits g — > and A — > and at the end the limit 
6 — s- 0. This procedure leads to the stability matrix 



B = 




(A.l) 



One of the corresponding eigenvectors is the unit 
vector in b direction and thus we are able to inves- 
tigate the behavior along this direction. The upper 
left block of this stability matrix is exactly the sta- 
bility matrix of the Einstein-Hilbert truncation and 
thus the solutions to the linearized flow equations 
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are known to be 

X k = (X ka -^ gka )(!f) 2 + ^ gka [l)\ 

g k = g k o(^) 2 - (A.2) 

Thus for the dimensionful cosmological constant A 
we find 

A fe - A fe0 + ^ (fc 4 - fc 4 ) . (A.3) 

Since A shows a k 4 scaling 5 of (3.12) is four as it 
was already derived in [37]. Thus we get the spectral 
dimension 

D s = £s = l- (A.4) 
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